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The exact amplitude for the asymptotic correlation function in the S=l/2 Heisenberg antiferro- 
magnetic chain is determined: 

< S$Sr >^ (-l) r 5 a6 (lnr) 1/2 /[(27r) 3/2 r]. 

The behaviour of the correlation functions for small xxz anisotropy and the form of finite-size 
corrections to the correlation function are also analysed. 
The asympototic behaviour of the equal-time correlation function in the Heisenberg antiferromagnetic chain has 
been difficult to determine numerically |[]-[5| because of the presence of a marginally irrelevant operator. This leads 
to a logarithmic factor of \/lnr and also to finite size effects which only vanish as 1/ In L where L is the system 
size. This marginal operator leads to logarithmic corrections, sometimes multiplicative and sometimes additive, to 
most long distance, low energy properties of the model. In particular recent experiments on Sr 2 Cu03 found evidence 
for the predicted [H logarithmic additive correction to the susceptibility ||. 

On the other hand, logarithmic corrections are absent for the xxz model, with Hamiltonian: 



H = + SfSf +1 + jStSf^}, 

i 

for 7 7^ 1. The model exhibits critical behaviour for — 1 < 7 < 1, with asympotic correlation functions: 



(i) 



G x (r) =< S£S* > 
G z (r) =< S^S Z > 



{-l) r A x r^ 
(-l)M.r- 1 /", 



with 



ri = 1 - [cos" 1 7]/tt (0 < i] < 1). 
What appears to be an exact formula for the amplitude A x (j) was recently conjectured Jul: 
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The main purpose of the present report is to determine the exact amplitude in the logarithmic, xxx case, 7=1, 
giving the result in the abstract. To do so it will be neccessary to consider the form of these correlation functions 
for 7 only slightly less than 1 where a crossover from logarithmic to non-logarithimic behaviour occurs. The other 
amplitude, A z , is not known in general. We will show that: 



lim A z /A x 

7->l- 



(6) 



The order of limits here is crucial; right at 7 = 1 the amplitudes of the (logarithmic) correlation functions G x and G z 
are equal. We also discuss the form of finite-size corrections for the correlation function on a ring of length L with 
periodic boundary conditions, G(r, L). 

The subsequent calculations are based on the continuum limit bosonized approximation to the xxz model. We 
follow the notation of H . The Hamiltonian density may be written: 



h = Ho (8n 2 /V3W(r L r R + jpi) + 9 z (j z L r R )}. 



(7) 



Here Ho is the Hamiltonian density for a free boson, of compactification radius R = l/\/2n, or equivalently, the SU(2) 
level 1 Wess-Zumino-Witten (WZW) non-linear a model. Jl,r are the left and right moving currents. (We set the 



1 



spin-wave velocity equal to 1.) For the isotropic model, with 7 = 1, g x — g z — g is of 0(1). The rather cumbersome 
normalization in Eq. (Q) is dictated by the convention that the operator multiplying g in the isotropic case have a 
correlation function with unit amplitude. For the xxz model with 7 close to 1, 

g*-g*<xl- r (8) 

These coupling constants obey the Kosterlitz-Thouless renormalization group (RG) equations: 

Z = dg z /d\nL = -(47r/V3)<& 

Px = dg x /d\nL = -(4ir/V3)g x g z . (9) 

The RG trajectories are sketched in Fig. |l[ g z (L) — g x (L) is an RG invariant along the flow. For g z > \g x \, the flow is 
to a fixed line, the positive g z axis. g1{L) — g x (L) — g1(oo) along these trajectories. Using the abelian bosonization 
formula J£ = — ( 1 />/8tt) (do + di)<p, we find that, at the fixed point, the effective Lagrangian is: 

C = (l/2)(d^) 2 [l - (27r S ,(oo)/V3]. (10) 
The staggered part of the local spin operators may be written in non-abelian bosonization notation as: 

S t a {-l) l trg<7, (11) 
where g is the two dimensional unitary matrix field of the WZW model. In terms of abelian bosonization: 

S*<x(-iy S in(ct>/R) 

Sf oc (-l) l cos(27ri?0), (12) 

where <j> denotes the dual field and R = \/\/2tt. From Eqs. (|l0|) and ([[2]) we can determine the correlation exponents 
of Eq. (|) with: 

77 = 1 - 2irg z {oo)/V3. (13) 

Note that, using Eq. (||), determined from the Bethe ansatz solution, the value of g z (oo) is determined exactly. 
The scaling dimensions of the staggered spin operators trga x and trga z are given by r//2 and l/2?7 respectively. To 
study the logarithmic behaviour, we will also need the anomalous dimensions for small non-zero gi, along the RG 
trajectories. These can be determined from the 3-point Green's functions < trga a J h L J^trga a > as in ^]. Using the 
fact that the operator product expansion gives: 

Jl{z)j\{z)g{z', z') - ^ )a \ b l2 +..-, (14) 
\Ztt(z — z ) y- 

we conclude that: 

< trger Q J b L J b R trga a >oc tr(a a a b a a a b ) = 4S ab - 2. (15) 
Thus, to linear order, the conclusion is: 

7x = l/2-(7r/V% z , 

7» = l/2+(7r/V3)(ff,-2 ffa! ). (16) 
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FIG. 1. The Kosterlitz-Thouless RG flows of Eq. (§). 



In discussing the asympototic correlation functions, using bosonization, it is convenient to introduce uniform and 
staggered terms: 

G\r) -> Gj,(r) + (-l) r G:(r), (17) 

where G^ and G* vary slowly on the scale of a lattice spacing. These two terms correspond to different Green's 
functions in the continuum limit field theory. In this paper we only discuss the staggered term. 
The staggered correlation functions (for an infinite spin chain) obey the RG equations: 

[d/dlnr + J2 PA9)d/d 9j + 2 7i (g)]G 4 s (r, g) = 0. (18) 

3 

Here g = (g z ,g x )- This follows from the fact that a rescaling of the length is equivalent to a change in the values of 
the effective coupling constants together with a rescaling of the fields with exponents 7 j. The solution of Eq. ([l8]) is: 

Gl(r,g°) = exp{-2 f dlnr , 7i 0(r , )]}F i 0(r)], (19) 

J r 

where the F l are arbitrary functions of g(r) 1 the solution of the RG equations, 

dg t /d\nr = f3 t {g). (20) 

Here g° = g(ro), denotes the value of the "bare" couplings at some reference "ultraviolet cut off" scale ro of order 
a lattice spacing. Since, for large r, g x (r) << 1, we may expand the functions Fi[g(r)] perturbatively in ^(r). 
Exactly this procedure is used to analyse deep inelastic scattering data in quantum chromodynamics. It is known as 
"renormalization group improved perturbation theory" . To lowest order these functions are just constants. Integrating 
the RG equations for the effective coupling constants, Eq. (||) we obtain: 

9x(i") = — — cosech(elnr) 

47T 

9z{r) = ^coth(elnr), (21) 

47T 

where we have defined: 

6 = 2(1-7?) =47r 3 ,M/V3. (22) 
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Now performing the integration over 7, (In r) in Eq. ( |l9| ) , we obtain: 



G x s (r) 
G z s (r) 



Ax (l - r -2eU/2 
r l-e/2 y X ' I 

A z (l-r-*y/ 2 



(23) 



r l+e/2 _|_ r -e)3/2 ' 

Note that we have defined the normalization constants so that the asympotic large-r behaviour is as in Eq. (g) 
Also note that, for e << 1, both correlation functions exhibit logarithimic behaviour over an intermediate range of r 
1 << lnr << 1/e. In this range of r, we obtain: 



G x 



G z s {r) 



(lnr) 



1/2 



V2eA x 

r 

V~eA z (lnr) 1 / 2 

2 3/2 r 



(24) 



Now consider taking the limit e — > 0, corresponding to the isotropic Heisenberg antiferromagnet. We see that in 
order for the correlation functions to remain finite at fixed r as e — > we must have A x oc 1/v^- Furthermore, in 
order to obtain the isotropic result, Gf (r) = Gf(r), we must have A z /A x — > 4, as e — > 0. Thus for small but finite 
e, Gf (r) Gg(r) in the intermediate range of r, but at very large r they exhibit slightly different exponents and 
amplitudes differing by a factor of 4. 

The exact amplitude, A x {rj) of Eq. (H) can be evaluated in closed form in the limit r\ — ► 1, e — ► 0. In this limit we 
may approximate sinh^t/ sinht « e _e */^ in the first term of the integrand and i] sa 1 in the second term. The integral 
can then be done exactly, giving: 



. 1 , 



1 



4 ( e ) 1/2^3/2- 

This diverges as 1/ ^/e, as expected. Thus we conclude, in the isotropic case: 

1 (lnr) 1 / 2 



G z s {r)=G*{r) 



(2tt) 3 / 2 



The asymptotic form of the Fourier transform for k w 7r is thus given by: 



G(A)= J] G W e ' 



3(2vr) 3 / 2 



ln|fc-7r|| 3 / 2 . 



(25) 



(26) 



(27) 



Note that the effect of the (lnr) 1 / 2 factor is to change the power of | In \k — n\\ from 1 to 3/2. If such a weak singularity 
could be observed, this formula might be useful to check the normalization in neutron scattering experiments. It follows 
from the above analysis that, for small xxz anisotropy, this isotropic formula remains valid down to exponentially 
small values of k — tt, making the log singularity of Eq. ( p7j ) observable. 

Several efforts have been made to check the field theory prediction of logarithmic behaviour numerically 
Hallberg et al. m obtained the above asympototic behaviour but with an amplitude of .06789 in place of the exact 
result (2tt)~ 3 / 2 = .06349364 . . .. This result was obtained from density matrix renormalization group calculations on 
rings of up to 70 sites using finite-size extrapolation. Koma and Mitzukoshi Q also obtained the above form with an 
amplitude of .065. [Alternatively, if they let the power of the logarithm be a free parameter they obtained a slightly 
better fit with a power of .47 instead of 1/2 and an amplitude of .071.] This was obtained using exact diagonalization 
results for L < 30 and zero temperature quantum Monte Carlo for 32 < L < 80. The agreement is remarkably 
good considering the severe difficulties of the extrapolation due to the logarithmic nature of the corrections. In the 
remainder of this report we consider the nature of the corrections to this formula, for the Heisenberg antiferromagnet. 

Let us begin with G s (r) for an infinite system. The integral in the exponent in Eq. ( |l9| ) can be rewritten as: 

/ h(g)/{3(9)}dg = (1/2) ln(r/r ) + / - + ]T a n g n 

** 90 J go |_ # n=0 

00 

= (l/2)In(r/r ) + (1/4) Hg(r)/g } + £ " 9o +1 }- (28) 



n=0 
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Here the a n terms arise from the higher order terms in the perturbative expansions of (3(g) and 7(g). Noting that 
all terms involving go ar e just constants, and also Taylor expanding the function F[g(r)] in Eq. (|l9|), we may finally 
write: 

1 00 

G fl (r)^— ==5> nff n (r), (29) 
r V9( r ) „=o 

in terms of some combined coefficients, b n . Including the cubic term in the /^-function for the isotropic case Jll]| : 

dg/d\nr = -(4tt/%/3) 5 2 - (1/2)(4tt/V3) 2 g 3 . (30) 



Integrating gives: 



Thus, we may write: 



-1- - 1 = (47r/V3){ln(r/r ) + (1/2) ln[ln(r/r )]} + 0(1). (31) 



rM 1 {ln(r/r ) + (1/2) lnHr/rp)]} 1 / 2 
s()= (2^r)372 r [l + 0(l/hir)]. (32) 



We may absorb the leading correction into a constant term inside the square root: 



°.w - „ ' { '° (Cr/r ° ) + (1/2)M '° Wro)1} ' /2 (i + o[i/(i»o']i. 

(27r) J /' i r 



(33) 



From Eq. (]31|) , C has the form: 

Q = £^3/4x90+0(1)^ (-34) 

The 0(1) term in the exponent in Eq. (|34|) could be computed. It requires calculation of the anomalous dimension 
j(g) to 0(g 2 ) and of the function F to 0(g). This term was ignored in Q leading to an inaccurate determination of 

So- 
Let us now consider the Green's function on a ring of length L, G s (r,r/L,g). The RG equation, Eq. fll8|), is 
still obeyed. The derivative in this equation may be taken either with respect to r or L with the ration r/L held 
fixed. This follows because a rescaling of both length scales is equivalent to a coupling constant redefinition. Using an 
L-derivative, the solution is now: 

G s (r,L, 5 °) = exp{-2 f dlnr' 7 [.g(r')]}^[ 5 (L), r/L], (35) 

The exponential factor is independent of r. The function F[g(L), r/L] may be expanded perturbatively in g(L) for 
large L: 

00 

F[g(L),r/L]=Y,9(LTFn(r/L). (36) 

The various functions F n (r/L) can be calculated by doing perturbation theory in the system with finite length. They 
should all obey the periodicity requirement: 

F n [r/L]=F n [(L-r)/L]. (37) 

If we take the asymptotic limit r/L — > 0, we should recover the infinite L result of Eq. (|3^). The zeroth order term, 
Fq(t / L) is obtained by ignoring the marginal interaction altogether and simply calculating: 

<tr(ag)(r)-tr(ag)(0) >L (38) 

in the conformally invariant WZW model, on a circle of length L. The correlation function on the circle (i.e. the 
cylinder in the space-time picture) is simply obtained by a conformal transformation and is given by: 
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< tv{ug){r) ■ tr(<? ff )(0) > L cx . (39) 

L sin(7rr/L) 



Thus we may write: 

G s (r,L) 



1 {ln(L/r ) + (l/2)ln[ln(L/r )]} 1 / 2 
(2tt) 3 / 2 (L/7r)sin(7rr/L) 



1 + Mk Awi » 



(40) 



for some other scaling function, F\ . Alternatively, solving the RG equation with an r-derivative, we obtaine this result 
with L replaced by r inside all logarithms and a different scaling function F[{r/L). [Note that, taking r » r$ with 
r/L held fixed, the difference between y/\n(r/ro) and yJ\n.{L/r ) is suppressed by a factor of l/rn(r/ro).] 

For the general xxz model the leading order finite-size scaling result is again obtained by the simple replacement: 

r — > (L/tt) sin(7rr/L). (41) 

In particular, for GJ in the xx model (7 = 0) we obtain: 

G x s oc [sin(7rr / 'L)]- 1 ' 2 . (42) 

The corrections are down by powers of 1/r rather than only logarithms. 

The efforts to fit numerical results on correlation functions in S=l/2 antiferromagets to a finite-size scaling form 
have a rather curious history. The case of Gf for the xx model was considered in Rather than using the result 
predicted by conformal invariance the authors adopted a phenomenological expression, with free parameters adjusted 
to obtain good data collapse, corresponding to the replacement: 



TTX 



l ' 2 



1 + . 28822 sinlf (1.673a;), (43) 

v sm TTx / 

for x = r/L. This leads to a correlation function not obeying the periodicity condition: 

G(r,L) = G(L-r,L). (44) 

Thus, the data fitting was only done for < x < 1/2. Over this range, these two functions actually agree to within 
about .05% as indicated in Figure |[ This indicates that the conformal field theory (CFT) prediction is extremely 
accurate for the xx model. It was proposed in (!]] that, in the general xxz model, one should use the form: 

[1 + .28822 sinh 2 (1.673a;)] 2 '', (45) 

for Gg. This is essentially the correct CFT prediction, due to the numerical agreement noted above. However, in jij] 
the exponent in Eq. (JSj) was taken to be a free parameter. For the Heisenberg model a best fit was obtained with the 
exponent 1.805 rather than the correct value of 2. Thus the scaling form used differed slightly from the one predicted 
by CFT as shown in Fig. § The maximum disagreement, at x = .5, is about 4%. 
Koma and Mizukoshi used the scaling function 

Gs(r ' X) ^ (LMMttt/L) ' (46) 

obtaining a best fit for A w .065 (close to (27r)~ 3 / 2 w .0635). The agreement between this formula and their numerical 
data is better than 1.26% for 1 < r < L/2 and 4 < L < 80. Taylor expanding in l/ln(L/ro), we see that this 
expression is consistent with Eq. ( fto| ) for a particular choice of the function F\, up to the small discrepancy in the 
amplitude. Eq. ( [l6"| ) has the great advantage of simultaneously having the correct periodicity property and the correct 
behaviour in the limit L — > 00. However, such an expression can only arise from Eq. ( |35| ) by summing an infinite 
number of terms in Eq. ( ]36| ) [and ignoring the ln[ln(L/ro)] terms in g(L)]. 

We expect that the somewhat larger discrepancy with CFT for the Heisenberg model than for the xx model can 
be accounted for by the log corrections. The range of r used in the numerical work of Hallberg et al. for which 
fairly good data collapse was obtained was only 10 < r < 30. In this range we might expect the factor l/ln(r/ro) in 
Eq. ( f4(i| ) (written with r replaced by L inside the logarithms) to be fairly constant. Thus the F[ term acts essentially 
as a small correction to the scaling function, ttx/ svn{irx). (A related observation was made in j^[.) It is feasible to 
push this renormalization group improved perturbation theory to one higher order and calculate Fi(r/L) in Eq. (|4C|). 
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This involves using the known result for the j3- function to 0(g 3 ), calculating the anomalous dimension to 0(g 2 ) and 
calculating the Green's function on a finite strip to O(g). We expect that this could give better agreement with the 
numerical results and could, in particular, reduce the small discrepancy between the exact amplitude and the results 
of 1 and §. 





FIG. 3. Comparison of the 2 different scaling functions for the xxx model. 
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